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ABSTRACT 



Context. Recent theoretical studies predict that the inner regions of galaxy clusters may have an enhanced helium 
abundance due to sedimentation over the cluster lifetime. If sedimentation is not suppressed (e.g., by tangled magnetic 
fields), this may significantly affect the cluster mass estimates. 

Aims. We use Chandra X-ray observations of eight relaxed galaxy clusters to investigate the upper limits to the effect 
of helium sedimentation on the measurement of cluster masses and the best-fit slopes of the Yx-Msoo and Yx-Afeoo 
scaling relations. 

Methods. We calculated gas mass and total mass in two limiting cases: a uniform, unenhanced abundance distribution 
and a radial distribution from numerical simulations of helium sedimentation on a timescale of 11 Gyrs. 
Results. The assumed helium sedimentation model, on average, produces a negligible increase in the gas mass inferred 
within large radii (r < rsoo) (1-3 ± 1.2 %) and a 10.2 ± 5.5 % mean decrease in the total mass inferred within r < rsoo- 
Significantly stronger effects in the gas mass (10.5 ± 0.8 %) and total mass (25.1 ± 1.1 %) are seen at small radii owing 
to a larger variance in helium abundance in the inner region, r < 0.1 rsoo ■ 

Conclusions. We find that the slope of the Yx — M500 scaling relation is not significantly affected by helium sedimentation. 

Key words. X-rays: clusters-galaxies: individual (MACS J0744.9+3927, MACS J1311.0-0311, MACS J1423.8+2404, 
MACS J1621. 6+3810, MACS J1720. 3+3536, Abell 1835, Abell 2204, Zwicky 3146) 



1. Introduction 



Clusters of galaxies are permeated by an optically thin ion- 
ized plasma at a temperature of ~ 10 8 K. This hot intra- 
cluster medium (ICM) is composed primarily of hydrogen 
and helium ions, with a small fraction (<1 %) of heav- 

■ ier elements, such as iron. A commonly adopted assump- 
tion is that helium and other heavy elements are uniformly 
distributed throughout the ICM, b u t theoret i cal st udies 

■ initiated bylFabian fc Pringld (|l977r i. iRephaelil (| 19781) and 
lAbramopoulos et al.l ()198lD suggest that diffusion processes 
will result in enrichment of iron and other heavy elements 
in th e cluster co r e. Fu r thermore, Gilfanoy fc Sunvaevl 
(fl98l. iQin fc Wul d2000l) . iChuzhov fc Nusserl (l2003l )aiid 
IChuzhov fc Loebl (12004TI have shown that the diffusion 
speed is greater for lighter elements, and conclude that 
the diffusion of helium ions into the core of galaxy clus- 
ters - helium sedimentation - can occur within a Hubble 
time. However the problem of helium sedimentation in 
galaxy clusters is still under debate. Recent theoreti- 
cal studies show that the presence of thermal diffusion 



Shtykovskiv fc Gilfanovl 120101 ) and strong magnetic fields 
Pen g fc Nagail 2009f ) in the intra-cluster medium may sup- 



press sedimentation of helium ions into the core of clusters. 



The diagnostics of abundance profiles in the intra- 
cluster plasma are based primarily on line emission from 
heavy elements. Since helium is fully ionized in clusters 
and produces no lines in the X-ray energy band, its radial 
distribution and abundance cannot be constrained by X- 
ray spectroscopic data alone. X-ray observations combined 
with observations of the Sunyaev-Zel'dovich effect (SZE) 
may be used to infer the helium a bundance in the intra- 
cluster medium (|Markevitchl 120071 ); however, the limited 
resolution and sensitivity of current SZE interferometers 
have so far not enabled the measurement of the distribu- 
tion of helium, even in the cluster cores where the sedi- 
mentation effect is the strongest. New capabilities for high 
angular resolution, high sensitivity SZE observations are 
currently being developed (e.g. at the Combined Array for 
Research in Millimeter-wave Astronomy, CARMA) which 
will allow observational constraints to be placed on helium 
sedimentation in cluster cores. 

The sedimentation of helium ions within clusters of 
galaxies may affect the physical quantities derived from 
X-ray observations such as cluster ga s mass, total mass 
(|Qin fc Wull2000l:lEttori fc Fabianll2006D . and t he cosmolog- 
ical parameters derived fro m these quantities (jMarkevitchl 
120071 : iPeng fc Na"gail l2009f ). Using high signal to noise 
Chandra observations of eight relaxed galaxy clusters, we 
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investigate a theoretical upper limit to the effect of helium 
sedimentation on the measurement of cluster properties; 
in particular, we study how this impacts mass determina- 
tions and the Yx-M^ scaling relations within overdensities 
A = 2500 and A = 500. To determine the upper limit of he- 
lium sedimentation effect we use a limiting case of a helium 
sedimentation profile for 11-Gyr-old galaxy clusters, based 
on th e numerical simulations provided by iPeng fc Nagail 
(2009). This work is the first application of a theoretically 
predicted helium sedimentation profile to Chandra obser- 
vations. 

This paper is organized as follows: in fJ3]we provide the 
details of our data reduction and modeling. The impact 
of helium sedimentation on X-ray derived gas and total 
mass and on the Yx-Ma scaling relation is presented and 
discussed in |21 In 2] we provide conclusions. In all calcu- 
lations we assume the cosmological parameters h — 0.73, 
Q M = 0.27 and n A = 0.73. 



Table 1. Cluster Sample 



Cluster 


z 


Nh a 
(cm" 2 ) 


Obs. ID 


Exp. Time 
(kscc) 


MACS J0744. 9+3927 


0.69 


5.6X10 19 


6111 


49.5 










3585 


19.7 










3197 


20.2 


Zwicky 3146 


0.29 


2.5 


xlO 20 


909 


46.0 










9371 


38.2 


MACS J1311. 0-0311 


0.49 


1.8 


xlO 20 


3258 


14.9 










6110 


63.2 










9381 


29.7 


Abell 1835 


0.25 


2.0 


xlO 20 


6880 


110.0 


MACS J1423. 8+2404 


0.54 


2.2 


xlO 20 


4195 


115.6 


MACS J1621. 6+3810 


0.46 


1.1 


xlO 20 


3254 


9.9 










6109 


37.5 










6172 


29.8 










9379 


29.2 










10785 


29.8 


Abell 2204 


0.15 


5.7 


X10 2Q 


7940 


72.9 


MACS J1720. 3+3536 


0.39 


3.5 


xlO 20 


3280 


20.8 










6107 


33.9 



(a) Lcidcn/Argcntinc/Bonn (LAB) Survey, see Kalbciia ct al. (2005) 



tions (unfiltered integration time > 60 ksec), all of which 
are available as Chandra archival data (see Table [lj. The 
sample spans a range of redshift (0.15 < z < 0.69). As 
part of the data reduction procedure, we apply afterglow, 
bad pixel and charge transfer inefficiency corrections to the 
Level 1 event files using C I AO 4 .1 and CALDB 4.1.1. We 
use the Markcvit ch et all (|2003h method for light curve 
filtering to eliminate flares in the background due to so- 
lar activity. The exposure times after filterin g are given in 
Table [J We follow the method described in Bulbul et al.l 
(|2010D for background subtraction. 



We use the 0.7-7.0 keV energy band extraction of spec- 
tra and images in order to minimize the effect of calibra- 
tion uncertainties at the lowest energies, and the effect of 
the high detector background at high energy. Spectra are 
extracted in concentric annuli surrounding the X-ray emis- 
sion centroid. In this process X-ray point sources were ex- 
cluded. An optically thin plasma emission model (APEC) is 
used, with temperature, abundance and normalization as 
free parameters in the spectroscopic fit (in XSPEC). The 
redshifts z and Galactic neutral hydrogen column densities 
Nh of the clusters are shown in Table [TJ We adopt a 1% 
systematic uncertainty in each bin of the surface brightness 
and a 10% system atic uncertainty in temperature profile 
(Bul bul et al.ll2010l) . The products of this reduction process 
are the X-ray surface brightness and temperature profiles 
shown in Figured) 




2. X-ray Observations and Modeling 

Our goal is to obtain density and temperature profiles for 
the gas in clusters based on Chandra X-ray imaging and 
spectroscopy, using a cluster model which accounts for a 
variable helium abundance as a function of radius. We de- 
scribe below the reduction of the X-ray data, the applica- 
tion of the cluster model, and the measurement of cluster 
gas mass and total mass. 



Fig. 1. The band-averaged X-ray cooling function, A ee (r), 
for two different helium distributions in the ICM: a uniform 
helium abundance (blue dashed line) and a helium sed imcn- 
tation distribution provided by IPeng fc Nagail ((2009) (red 
line) . The overdensity radius r5oo is the radius within which 
the mean cluster density is a factor of 500 times greater 
than the critical density of the universe at the cluster red- 
shift. 



2.1. X-ray Data Reduction 

We ch ose the eight galaxy clusters from the lAllen et aLI 
(2008) sample with the deepest Chandra ACIS-I observa- 
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2.2. Application of Cluster Models to X-ray Data 

We start with the projected surface brig htness 
(Birki nshaw et al.lll99ll : iHughes fc Birkinshawl[T998h : 



S x — 



1 



47r(l + zf 



n\ A ee d£, 



(1) 



where S x is in detector units (counts s 
A ee is the band averaged X-ray cooling function in counts 
cm 3 s , n e is the number density of electrons, z is the 
redshift, and the integral is along the line of sight, £. 

We calculate the X-ray cooling f unction in two lim- 
iting cases: (i) a uniform abundance (lAnders fe Grevesse 
1989) and (ii) the radial abundance profile of lPeng fe Naeai 
(2009) obtained for a sedimentation time scale of 
11 Gyr with no magnetic field. We use the astro- 
physical plasma emission database (APED) which in- 
cludes thermal bremsstrahlung, radiat ive recombination , 
line emission, and two-photon emission (|Smith et al.ll200ll) . 
In addition, we apply correc t ions for galactic absorp- 



tion (iMor 



(Itoh ct al 



Itoh et al 



ison fe McCammonl Il983f ). relativistic effects 
2000), and electron-electron bremsstrahlung 



2001), to obtain the X-ray emissivity (e v ) as 
a function of frequency v. 

The band-averaged X-ray cooling function (A ee (T, r)) is 
obtained by integrating the e„ over the observed Chandra 
energy band (0.7-7.0 keV): 



Ke{T,r) = 



(2) 



The resulting band- averaged X-ray cooling fu nction for a 
uniform abundance (lAnders fe Grevesse! 19891) and a he- 
lium sedimentation profile dPeng fe Nagail l2009h is shown 
in Figure [T] 

To obtain the electron number density n e (r) and tem- 
perature T{r) profiles, we apply the analytic model for 
the intra-cluster plasma developed by iBulbul et al.l ([20101 ) 
and compare the model predictions to the observed surface 
brightness and temperature data points found in H2.ll 

The electron number density is 



n e (r) =n 



r() 



(l + r/r s f- 



08-2) r/r s (l + r/r s y- 



;» ( 3 ) 



where n e Q is the central value of the number density of 
electrons, r s is the scaling radius, (3 + 1 is the slope of the 
total matter density profile, and n is the polytropic index. 
Tcooiir) is the phenomenol ogical core-tape r function used 
for cool-core clusters ( Vikhlinin et al.|[2006l ). 



i(r) 



a + (r/r coo/ ) 7 



(4) 



1 + (r/r coo/ )T ' 

The corresponding lBulbul et all (|2010l ) temperature profile 



is 



T(r) =T 



1 



(l + r/r s ) 



P-2 



1 



08 -2)r/r s (l + r/r s y- 2 



oi(r) (5) 



where To is the normalization parameter, and the other 
parameters are in common with those in Equation [31 

We vary the model para meters using a Monte 
Carlo Markov Ch ain approach (|Bonamente et af] [2004; 
Bulbul ct al.l |2010T) ; the best fit model parameters for a 



uniform helium abundance and a helium sedimentation sce- 
nario are given in Table [3] and are shown by the solid blue 
and red curves in Figure [3J We assume that clusters with 
redshift z > 0.3 have /3 = 2, since the polytropic index 
(n) and j3 cannot be determined simultaneously from the 
data available. For all of the clusters we use a core taper 
parameter 7 — 2.0 except for Zwicky 3146, for which we 
obtained a better fit with 7 = 1.0. With these choices of 
fixed parameters, we obtain acceptable fits to the data for 
both the assumed uniform and sedimented helium profiles 
(see Tables Hand O. 

2.3. Cluster Mass 

The gas mass is the volumetric integral of the number den- 
sity of electrons multiplied by the electron mean molecular 
weight: 



M gas (r) — Trip I ^i e (r) n e (r) dV, 
Jv 



(6) 



where m p is the proton mass and n e (r) is defined in 
Equation [3] Since the mean molecular weight of electrons 
/Lt e (r) is dependent on the ion distribution in the plasma, 
we use the same approach as in the X-ray cooling func- 
tion calculations to determine its radial distribution. The 
electron mean molecular weight n e jr) , which is calculated 
for a uniform helium abundance (|Anders fc Grevessdll989f) 



and t he radial helium abundance profile of (jPeng fc Nagail 
2009) obtained for a sedimentation timescale of 11 Gyr with 
no magnetic field, is shown in the top panel of Figure [3J 
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o 
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Fig. 3. The distributions of the electron mean molec- 
ular weight (top panel) and the total mean molecular 
weig ht (bottom panel) a s a fun ction of radius for a uni- 
form lAnders fc Grevesse] (|1989T) helium abundance (blue 
dashed line) and for t he 11 Gyr helium sedimentation model 
iPeng fc Nagail (|2009| )(red line). 
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Fig. 2. X-ray surface brightness (left column) and temperature profiles (right column) for the sample of eight relaxed 
galaxy clusters. The black points are derived from the X-ray i mage and spectroscopic dat a. The blue lines in both profiles 
show the best fit model to the data obtained us ing a uniform Anders fc Grevessel (|1989|) helium abundance, and the red 
lines show the best fit model obtained using the iPeng fc Naeai ( 2009t ) helium sedimentation distribution; the green lines 
in surface brightness profiles indicate the background levels determined from the blank sky observations. The overall x 2 
of the fits obtained from the uniform helium abundance and the helium sedimentation model (Tables [5] and [3]) show that 
Doth distributions give equally good fits to Chandra data. 
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Similarly, the total mass is the volumetric integral of 
the total matter density /3t t(r), 



(7) 



M tot (r) = / p tot (r) dV. 



As in I Bulbul et al.l (|2010() . the total matter density is found 
by 



Ptot(r) 



T fc(n + l)(/3-l) 



Htot(r) (1 + r/r s )r/r s 



where k is the Boltzmann constant, G is Newton's gravi- 
tational constant, and /it t(r) is the total mean molecular 
weight, which is also dependent on the ion distribution in 
the plasma. The total mean molecular weight, [itot(r) is 
calculated assuming a uniform helium distribution and the 
helium sedimentation model (see Figure [3] bottom panel) 
as was done for the [i e calculations. 



3. Results 

We measure cluster masses and scaling relations as- 
suming the limiting case s of uniform helium abundance 
(| Anders fc GrevesseTll989D and the radially- varying sedi - 
mented helium abundance profile of iPeng fc Nagai (|2009T ) 
described in Section! 



In order to determine the fractional change in cluster 
mass measurements introduced by the spatial variation in 
the helium distribution, we calculate gas mass and total 
mass for each helium abundance model. Gas mass is cal- 
culated using Equation H2 and the electron mean molecular 
weight distributions shown in the top panel of Figure |3] 
The total mass is calculated using Equation [7] and the total 
mean molecular weight distributions shown in the bottom 
panel of Figure [3] The gas mass and total mass at r25oo 
and rsoo are given in Tables H] and [5] for both the uniform 
helium abundance and the sedimented case. 

The fractional change with radius in gas mass and to- 
tal mass measurements as a result of an extreme helium 
sedimentation case is shown in Figure 0] for each of the re- 
laxed galaxy clusters in our sample. The black line shows 
the mean fractional change in gas mass (left panel) and 
total mass (right panel). The shaded areas indicate the 
standard deviations of the gas mass and total mass in- 
ferred for the sample. In Table |6] we show the weighted 
mean percentage difference with rms errors in gas mass 
and total mass measurements within small cluster radius 
(r < O.lrsoo), intermediate radius (r < T25oo) an d large 
radius (r < r 500 ). On average, at large radii (r < r^oo) 
the effect of helium sedimentation on gas mass is negligible 
(1.3 ± 1.2 %). The helium sedimentation model produces 
an average of 10.2 ± 5.5 % decrease in the inferred total 
mass within r$oo- Significantly stronger effects in the gas 
mass (10.5 ± 0.8 %) and total mass (25.1 ± 1.1 %) are 
seen at small cluster radii (r < O.lrsoo) where the helium 
abundance enhancement and gradient are greater. This re- 
sult is in agreement with the expected change i n the mass 
measurements predicted b y theoretical studies (|Qin fc Wul 
120001: IPeng fc Nagaill2009h . 



3.1. Y x -M Scaling Relation 

We examine the effect of helium sedimentation on the clus- 
ter mass s caling relations. The X-ray mass proxy, Yx, is 
defined as (jKravtsov et al.ll2006l ). 



Y : 



x 



M gas (r 500 ) T X - 



(9) 



where Tx is the me an X-ray spectroscopic temperature 
meas ured within rsoo (|Kravtsov et al . 2006; Vikhl inin et al.l 
2009). The importance of this quantity is that there is 
a corr elation between Yx a nd the total cluster mass, 
M 500 (jKravtsov et al.l 120061: IVikhlinin etHl I2009D and 
this correlation motivates Yx as a proxy fo r total mass 
()Vikhlinin et al.l 120091: lAndersson et al.ll2010f ). In this sec- 
tion we investigate the effect of helium sedimentation on the 
Yx-Ma scaling relation at overdensities A = (2500, 500) 
(i.e. for properties computed within r < (^2500, r 50o))- 

The Yx-Msoo scaling relation is found using gas mass 
and total mass measurements shown in Tables 2] and [5] 
for a uniform helium abundance and helium sedimenta- 
tion distributions. The average spectral temperature Tx is 
obtained with a single-temperature thermal plasma model 
(APEC) fit to the Chandra spectrum in the radial range 
0.15r 50 o < r < r 5no . Using the XSPEC analysis package 
to vary the helium abundance assumed in our spectroscopic 
fits, we determine that the assumed helium distribution has 
a negligible effect (< 1%) on the measurements of the av- 
erage X-ray spectral temperatures Tx- We report the aver- 
age X-ray temperature Tx and the X-ray mass proxy Yx 
for a uniform helium distribution and helium sedimentation 
model in Table [7J 

The Yx-M 50 o scaling relation for our sample of eight 
relaxed galaxy clusters is shown in Figure [3] (left panel). For 
each A, we fit the Yx -M& data using a power law relation 
(|Vikhlinin et al.ll2009ft . 



M A E{z) 



2/5 



x 



3 x l0 li M Q keV 



h 1,2 M t 



01 



(10) 



where a is the normalization, b is the slope, and E(z) is the 
evolution function. The best fit normalizations and slopes 
are reported in Table |8] with the goodness of the fit. In 
Figure[5]the dashed lines shows the best fit power law model 
and solid lines shows the 90 % confidence intervals for both 
a uniform helium abundance and the helium sedimentation 
cases. 

From the Yx-M 500 fit to the data on the eight galaxy 
clusters in our sample, we conclude that both the uniform 
and sedimented helium fits are statistically con s istent with 
the power law slope found by IVikhlinin et al.l (|2009f ) (see 
Table [5]) . helium sedimentation changes the normalization 
by 10 ± 0.9 % at r 50 o- The difference in best-fit slope 
between the uniform and helium sedimentation models is 
small compared to the statistical uncertainty. 

We also investigate the effect of helium sedimentation 
on the Yx-M 2 5oo scaling relation. Figure [5] (right panel) 
shows the Yx-M 2 5oo scaling relation for our sample. We fit 
the YX-M2500 scaling relation data with a power law shown 
in Equation [TU] The dashed lines show the best fit power 
law relation with the goodness of fit reported in Table [9] for 
the uniform and sedimented helium profiles. In Figure [3] the 
best fit power law model is shown in dashed line and 90% 
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Fig. 4. The fractional change in gas mass an d total mass with radius as a result of an extreme helium sedimentation 
with respect to masses produced by a uniform [Anders fc Grevessel (|l989l ) helium abundance for a sample of eight relaxed 
clusters. The black lines correspond to the mean fractional change and the shaded are indicates the standard deviation 
of the fractional change in gas mass and total mass measurements. 
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Fig. 5. The Yjf-Msoo and Y^-M^oo scaling relations for a sample of eight relaxed galaxy clusters. The blue data points 
and curves correspond to the measurements obtained from a uniform helium abundance case and red data points and 
curves correspond to the measurements obtained from a helium sedimentation model. The dashed lines show the best fit 
power law relation and solid curves show the 90% confidence levels. 



confidence levels are shown in solid lines for both uniform 
helium abundance and helium sedimentation models. 

The self-similar Yx — M 2500 scaling relation discussed 
by iBonamente et all ( 200811 is of the form: 



Y X oc M^ E(z) 2 ^ 



which is equivalent to the Yx 
form we use for this work 



(11) 

-M2500 scaling relation 



M 25m E{zf' b «yf . 



(12) 



Recently, IBonamente et al.l (|2008l) measured the Yx — 
M2500 scaling relation, and found it to be consistent with 
the self-similar prediction, b = 1.66 ± 0.20 in Equation [TOT 
In our fit Yx — M2500 (shown in Table , we also find that 



both the uniform helium abundance and the helium sedi- 
mentation cases are con sistent with the slope measured by 
IBonamente et all (|2008l) and with the self-similar expecta- 
tion (Equation IT2"j) . The difference in best- fit slopes b is 
negligible compared to the statistical uncertainty reported 
in Table [9l However, we find that helium sedimentation af- 
fects the normalization of the Yx — -M2500 scaling relation 
by 13 ± 0.6 %. 



4. Conclusion 

In this paper we investigate the upper limits to the ef- 
fect of the helium sedimentation on X-ray derived cluster 
masses and the Yx — M500 and Yx — M2500 scaling rela- 
tions using Chandra X-ray observations of eight relaxed 
galaxy clusters. We used a limiting helium sedimentation 
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profile for 11 Gyr old clus t ers ba sed on the simulations per- 
formed bv lPeng k Naerail (|2009t ) and compare these results 
with those assuming a uniform helium distribution in order 
to determine the maximum impact of helium sedimenta- 
tion. This work is the first application of a theoretically 
predicted helium sedimentation profile to Chandra obser- 
vations. We fit the deep exposure Chandra X-ray spectro- 
scopic and imaging data of the eight relaxed galaxy clusters 
with an analytic model of the intra-cluster plasma devel- 
oped bv lBulbul et al.l (|2010h . We demonstrated that both a 
uniform helium and a limiting helium sedimentation model 
can accurately describe the surface brightness and temper- 
ature profiles obtained from the Chandra X-ray data. 

We have found that, on average, the effect of helium 
sedimentation on gas mass inferred within large radii (r < 
rsoo) is negligible (1.3 ± 1.2 %). The helium sedimentation 
model produces a 10.2 ± 5.5 % mean decrease in the total 
mass inferred within r < r$oo. Significantly stronger effects 
in the gas mass (10.5 ± 0.8 %) and total mass (25.1 ± 1.1 
%) are seen at small radii owing to a larger variance in 
helium abundance in the inner region, r < O.lrsoo- This 
study supports the view that helium sedimentation should 
have a negligible impact on cluster mass inferred within 
large radii. This result is consisten t with the predict i ons o f 
the previous theoretical s tudies by lEttori fc Fabianl (l2006h 
and lPeng k Nagail (|2009fl . 

The fractional change in both the gas mass and the 
total mass measurements due to helium sedimentation do 
not show any trend with redshift. The strongest effect on 
the cluster total mass is observed on an intermediate red- 
shift cluster, Zwicky 3146, while the smallest effect is ob- 
served on the highest redshift cluster in the sample, MACS 
J0744.9+3827. 

In order to investigate the effect of helium sedimen- 
tation on the best-fit slopes and normalizations of the 
Yx — M 5no and Yx — M2500 scaling relations, we used 
gas mass measurements, and the mean X-ray spectroscopic 
temperature Tx to estimate the X-ray mass proxy Yx 
for both the uniform and sedimented helium abundance 
profiles. Both uniform and helium sedimentation models 
produce slopes which are st atistically consistent wit h the 
power law slop e s fou nd by iBonamente et al.l ((2008) and 
Vikhlini n et al.l ((2009:). We have found that helium sedi- 
mentation has a negligible effect on the slopes of Yx — -M500 
and Yx — M2500 scaling relations. We have also found 
that helium sedimentation changes the normalization of the 
Yx— M500 scaling relation by 10 ± 0.9 % and the Yx— Af 2 5oo 
scaling relation by 13 ± 0.6 %. 
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Table 2. Best- fit Parame ters of the iBulbul et al.l (|2010T) Model Produced Using a Uniform helium Abundance Model 
(|Anders fc Grevessdll989t) 



Cluster 


n c0 
(10~ 2 cm~ 3 ) 


r s 
(arcsec) 




n 







To 
(keV) 


(arcsec) 




CK 




7 


X z (d.o.f.) 


P value 


MACS J0744. 9+3927 


rj i q+0.69 

^' ici -0.30 


in ^ * +2. 14 

13.64+.3 18 


2. 


__+0. 15 
' '-0.18 


2.0 




22.47±?-™ 


QQ + 4.56 


0. 


i i+0 

!3+o 


08 
.02 


1.0 


44.3 (51) 


73.5% 


Zwicky 3146 


i n— +0 43 


i-r-i nn+5.47 

52 -"I 5 .43 


2 


ro + 0.83 
■ 58 Io.77 


2.85_q 


83 
.18 


30 - 02 -2.31 


i-n __r--t-5. 99 

58.95+°™ 


0. 


06 -0 


02 
.02 


1.0 


66.8 (78) 


81.2% 


MACS J1311. 0-0311 


, -O + 0.20 
J- 1 ° — 0.65 


40.32±|;|» 


5. 


oc + 3.19 
oo _0.47 


2.0 




U.36±S;« 


20.73l 4 j° 


0. 


33+° 


04 
.13 


2.0 


70.56 (63) 


42.5 % 


Aboil 1835 


„ i. 7 +0.29 
z '°'-0.07 


AO qh + 3-29 
W - OZ -6.48 


3 


qo + 0.73 


1.94t° 


15 

.22 


18.26±?;£ 


22.65+_5:? 7 


0. 


18±° 


02 
.01 


2.0 


99.3 (93) 


30.8% 


MACS J1423. 8+2404 


, 9Q + 0.68 
°- zy -0.29 




2 


yo -0.16 


2.0 




16.03±i;~ 


13.07l° ; 7 3 


0. 


18±° 


04 

.02 


2.0 


45.5 (50) 


65.3% 


MACS J1621. 6+3810 


2 84+ 62 

^•°*-0.16 


14.39l 3 2 ;f 5 


3 


10 +0.17 
0.12 


2.0 




ia.«»±S:S? 


9 oi+°' 73 


0. 


26l° 


05 
.03 


2.0 


32.4 (45) 


92.1% 


Abell 2204 


^ z -0.24 


21.73^1° 


6. 


44 to.5? 


1.39±° 


04 
.06 


14.28±° ; 7 | 


19 42+ 60 


0. 


16l° 


01 
.01 


2.0 


115.5 (145) 


96.6% 


MACS J1720. 3+3536 


„ ,„ + 0.24 
Z '- L -0.20 


31.69±1:S1 


3 


Q7 +0.29 
y ' -0.27 


2.0 




13.011°- 


10 95+ 86 


0. 


13-°o 


02 
.01 


2.0 


46.5 (60) 


89.9% 



Tabl e 3. Best-fit Parameters of the lBulbul et al.l (|2010f ) Produced Using a helium Sedimentation Model (jPeng fc Naeail 
2009) 



Cluster 


(10~ 2 cm -3 ) 


r s 

(arcsec) 




n 


P 


To 
(keV) 


(arcsec) 




a 


7 


X 2 (d.o.f.) 


P value 


MACS J0744. 9+3927 


2 07+ - 26 

z - u ' -0.25 


14.79til 7 


2 


g4 +o.io 

°^-0.15 


2.0 


21.751;- 


1 2 1 q+ 3 - 51 





.13 


1.0 


45.7 (52) 


71.85% 


Zwicky 3146 


, , 8 +0.10 
1 - 2ii -0.07 




2 


fi o + 0.39 
bS -0.24 


9 ca+0.37 
2bS -0.13 


25.92l ; 33 


47 33+ 3 63 

"• M -5 .62 





.11 


1.0 


60.83 (79) 


93.6% 


MACS J1311. 0-0311 


J-.OD_0.09 


42.26± 6 7 ; 4 6 


5 


OO + 0.54 
oo -0.67 


2.0 


10 91 + 101 


19.3211;?* 





.34 


2.0 


44.6 (64) 


68.9% 


Abell 1835 


9 oc+0.30 
-'■OO-O.OS 


38.76+«;}° 


4 


n1 +0.28 
,u — 0.83 


o oo"t"0' 35 
2-88_ 29 


17.75i° ;« 


21.43+°% 





.18 


2.0 


108.1 (94) 


13.6 % 


MACS J1423. 8+2404 


2 99+ 11 


, , -, + 0.86 
J. 1 ' 'O-0.84 


2 


03+0.11 

MJ -0.54 


2.0 


I5.49i;;» 


12.28l° .S 





.18 


2.0 


47.2 (51) 


58.5 % 


MACS J1621. 6+3810 


2 ro+0.16 
z - oo _0.14 


14.61+i-ig 


3 


g + 0.09 
UD -0.12 


2.0 


12.661?;?] 


8.29+g;~ 





.26 


2.0 


33.88 (46) 


88.8 % 


Abell 2204 


4 15+ - 07 

^■J-O_ .08 


19.8X+S-.JS 


5 


CQ+0.66 
-0.39 


2 4 q+o.04 

z -^°-0.05 


14 51+ - 53 

J^-OJ_0.59 


18.81±g;»S 





.16 


2.0 


127.7 (146) 


84.6 % 


MACS J1720. 3+3536 


1 91+ 013 

J- y -0.12 


31.62t 3 5 r 3 l 


3 


85+ 21 


2.0 


12 96+ 101 
i_.ao_ 79 


10.57l° 6 6 6 3 





.13 


2.0 


47.3(61) 


90.0% 
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Table 4. Gas Mass and Total Mass for the Uniform helium Abundance Model at the Overdensity Radii r 2 5oo and r 50 o 



Cluster 


7^2500 

(arcsec) 


(1O 13 M ) 


tot \ ' Z OUU / 

(1O 14 M ) 


/"500 

(arcsec) 


Mnn *(r<\()()) 

(10 13 M Q ) 


(1O 14 M ) 


MACS J0744.9+3927 


57.9±li 


2.98±°-JS 


9 r»7+0.36 
z - u ' -0.29 


120.2± 8 6 j 


8.78tS;$j 


fi Q+0.80 


Zwicky 3146 


127.6±S:| 


4 26+ ' 18 

^• ZD -0.20 


o i 7 +0.35 
-0.38 


266.7t 1 4.2 


n 79+O.39 
y - ' z -0.44 


5.781!;^ 


MACS J1311.0-0311 


78.7+j-j 


2 29+°' 19 

z - za -0.18 


o ci +0.61 
Z - O1 -0.49 


173.1+-;! 


r 1fi +0.29 
d - io -0.29 


5-34±!i| 


Abell 1835 


150.61^2 


4: - al -0.17 


o 7 9+0.26 
°- ' -0.30 


sog^t^ 8 ! 


12.08+S; 38 


6.38tr 7 9 


MACS J1423. 8+2404 






1 61+ ' 19 


125.9±g;l 


5.46i°j 3 7 


9 C9+0.34 

Z - Oz -0.37 


MACS J1621.6+3810 


67.8±|-J 


i-78±8:i| 


i 07+O.27 
1 - ot -0.20 


135.1±?;1 


4.77t°; 3 6 


9 1 7+0.48 
Z - Ll -0.35 


Abell 2204 


225.711! 


o qq+0 09 


q O7+O.I9 


479.8+!!'j 


10.35t°; 2 J 


6.48±r 4 I 


MACS J1720.3+3536 


9i-6±|:l 


9 09+O.I8 
z - oz -0.16 


9 0,4+0.34 
^•O1_o.28 


190.5± 9 ° 3 6 


7" QQ+0.35 
< - oo -0.33 





Table 5. Gas Mass and Total Mass for the helium Sedimentation Model at the Overdensity Radii r 250 o and rsoo 



Cluster r 2 5oo M 9as (r 25 oo) M to t(r 25 oo) r 5 oo M gaa (r 500 ) M tot (r 5 oo) 
(arcsec) (1Q 13 M ) (10 14 M Q ) (arcsec) (1O 13 M ) (10 14 M Q ) 



MACS J0744 


.9+3927 


53.9+J.I 


9 7C+O.21 
z - ' J -0.22 


1. 


fifi +0.27 
• DD -0.27 


116.9+1° 


8.43+°,;*° 


4 i+0.56 
■J- 4i -0.56 


Zwicky 3146 




116.0+ 4 ; 2 . 


04+0.17 

°- M ^-0.15 


2 


OO+0.27 
•<*°-0.22 


246.5±!?; 1 , 


9.14l«; 44 


4 56+ - 77 

^•OO_0.61 


MACS J1311 


.0-0311 


71.2+ 3 j 


9 19+O.I2 
z - ±z -0.15 


1. 


o fi +0.32 
° D -0.32 


165.6+? 2 5 


5.06t°;l| 


4 fi 7+0.82 
-0.83 


Abell 1835 






4 62+ ' 10 

4 - DZ -0.17 


2 


oq+0.16 


294.8+! 2 7 7 


11.59+ ,; 3 ? 


r c-o+0.45 
J-JO-0.69 


MACS J1423 


.8+2404 


58.8 + ^2 


2 01+ 010 
z - ul -o.oi 


1. 


9O+0.17 
• zy -0.14 


120.9+. 4 ; 8 , 


5 20+ ' 19 

d - zu -0.18 


9 90+O.27 
Z - ZO -0.23 


MACS J1621 


.6+3810 


63.41 3 ; 4 


i - DO -0.11 


1 


19+O.19 
• iz -0.17 


130.7+6;! 


4-55i«; 23 


1 Q7+0.33 
l-y_0.29 


Abell 2204 




201.0t 3 j 


cc+0.08 
O.OO_ 09 


2 


OO+0.13 
• J0 -0.14 


440.4t! ; 5 9 


9.53±° 2 £ 


5 01 + ' 40 

°- ul -0.36 


MACS J1720 


.3+3536 


84.6± 3 ;?. 


2.59«;!« 


1 


04+O.27 
• 84 -0.22 


182.9lg;i 


-09 


79+O.54 
°- ' z -0.47 



Table 6. Percentage Changes in M gas and M to t 

Region % Change in M gas % Change in M tot 

0.0 < r < 0.1r 500 10.5 ± 0.8 25.1 ± 1.1 

0.0 < r < r 250 o 3.5 ± 1.0 13.9 ± 3.1 

0.0 < r < r 500 1.4 ± 1.2 12.5 ± 5.5 
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Table 7. Measured Average Temperature (Tx) and the Mass Proxy (Yx) 



Cluster 


T x 
(keV) 


Y x 

(Uniform heliumAbundance) 
(lO 14 fceVM ) 


Y x 

( helium Sedimentation) 
(10 14 fceVM Q ) 


MACS J0744.9+3927 


6.80 ± 0.5 


5.97 ± 0.58 


5.73 ± 0.54 


Zwicky 3146 


7.29 ± 0.4 


7.09 ± 0.48 


6.66 ± 0.45 


MACS J1311. 0-0311 


5.75 ± 0.7 


2.97 ± 0.39 


2.91 ± 0.37 


Abell 1835 


8.79 ± 0.3 


10.57 ± 0.49 


10.30 ± 0.50 


MACS J1423.8+2404 


6.11 ± 0.5 


3.34 ± 0.31 


3.18 ± 0.27 


MACS J1621. 6+3810 


6.93 ± 0.5 


3.31 ± 0.29 


3.15 ± 0.27 


Abell 2204 


6.91 ± 0.3 


7.16 ± 0.36 


6.63 ± 0.34 


MACS J1720.3+3536 


7.63 ± 0.8 


5.63 ± 0.64 


5.41 ± 0.61 



Table 8. Comparison of Normalization and Slope of the Yx — M 500 Scaling Relation 

a b x z (d.o.f) F(%) 

(10 14 h 1/2 Mq) 

Uniform helium Abundance 3.12 ± 0.87 0.83 ± 0.42 8.6 (6) 19.7 
helium Sedimentation Model 2.92± 0.61 0.77 ± 0.22 10.7 (6) 9.8 



Table 9. Comparison of Normalization and Slope of the Yx — M2500 Scaling Relation 





a 

(10 14 h 1/2 Mq) 


b 


X' (d.o.f.) 


P 


Uniform helium Abundance 


1.94 ± 0.50 


0.72 ± 0.28 


5.0 (6) 


54.3 


helium Sedimentation Model 


1.66 ± 0.30 


0.66 ± 0.11 


5.6 (6) 


46.9 



